ISME Communications

The Deep Chlorophyll Maximum (DCM) is likely the most important feature organizing the marine epipelagic environment. Within this layer, opposing gradients of light and nutrients create a stratified habitat that supports high phytoplankton biomass and a substantial fraction of oceanic primary production. Despite its ecological importance, most studies treat the DCM as a single depth, overlooking its fine-scale heterogeneity. Here we investigated prokaryotic community organization across the DCM in the northwestern Mediterranean Sea through high-resolution sampling of four profiles collected over two days. Free-living (0.2-3 {micro}m) and particle-associated (3-20 {micro}m) communities were characterized using 16S rRNA gene amplicon sequencing. Prokaryotic communities changed progressively along the vertical gradient, revealing the DCM as a microbial coenocline with continuous community turnover. Fuzzy clustering identified distinct assemblages associated with environmental transitions from warm surface waters to the chlorophyll maximum, the nitrite peak below the DCM, and deeper nitrate-rich layers. In both the free-living and particle-associated fractions, most ASVs remained consistently associated with the same depth-defined clusters across all samplings, indicating stable niche partitioning over short timescales. However, these temporally stable ASVs accounted for a substantially smaller fraction of community sequences in particle-associated communities, suggesting higher dynamism, likely driven by particle-mediated transport. Nevertheless, phylogenetic analyses revealed that closely related ASVs tended to occupy similar depth niches, indicating that habitat preferences are phylogenetically conserved in both size fractions. Our results demonstrate prokaryotic niche partitioning over scales of only a few meters within the DCM, highlighting the importance of fine-scale sampling for understanding microbial community structure and responses to ocean change.

Sialic acids (Sias) and related nonulosonic acids are critical components of glycoconjugates involved in host-pathogen interactions, immune regulation, and cell signalling. Despite their biotechnological relevance, the diversity of enzymes involved in Sia biosynthesis remains largely underexplored due to limitations in culture-dependent methods and the lack of (ultra)high-throughput screening strategies. Here, we report the development of a highly sensitive droplet-based microfluidic screening platform enabling the functional discovery of sialic acid aldolases in environmental metagenomes. The method integrates a fluorescence-coupled enzymatic cascade compatible with fluorescence-activated droplet sorting (FADS), allowing the screening of >10 droplets per experiment, as well as a downstream validation strategy for the selected hits. Although some limitations were identified, the system demonstrated high sensitivity and was utilised for the screening of a metagenomic library from garden soil. During this campaign, a potential new sialic acid aldolase enzyme was identified. This work establishes a generalizable framework for measuring complex, multi-step enzymatic functions at ultrahigh throughput using coupled cascades in droplets

Aerobiological studies have largely focused on near-surface sampling and horizontal biogeographic patterns, while vertical structuring of airborne microbial communities within the atmospheric boundary layer (ABL) remains poorly understood. Here, we investigated microbial communities across the lower and upper ABL in a low-orography coastal site on the Antarctic Peninsula, representative of the Southern Ocean marine ABL and with low direct human influence. Airborne microorganisms were sampled simultaneously using ground-based and aerial platforms on five occasions. Community composition, abundance, and cell morphometry were assessed using metabarcoding and epifluorescence microscopy and interpreted alongside atmospheric observations. Airborne bacterial and eukaryotic communities showed consistent vertical stratification, although partial taxonomic overlap indicates vertical connectivity between atmospheric layers. Lower ABL communities were more diverse than upper ABL counterpart, compositionally homogeneous, and dominated by marine-associated taxa, reflecting strong influence from local sources and turbulent mixing. In contrast, upper ABL communities were less diverse but more heterogeneous among sampling events, enriched in stress-tolerant, terrestrial and plant-associated taxa, consistent with atmospheric filtering, selective upward transport, and long-range atmospheric inputs. Upper-layer samples also exhibited higher microbial abundance and greater prevalence of elongated cell morphologies, suggesting particle accumulation aloft and aerodynamic selection permanence. Together, these findings identify the Southern Ocean ABL as a vertically structured microbial habitat organized into two partially decoupled sublayers, in which atmospheric dynamics regulate microbial dispersal, ecosystem connectivity, and biogeographic patterns.

Phage tail-like elements (PTEs) -- tailocins, bacterial type VI secretion systems (T6SS), and extracellular contractile injection systems (eCIS) -- are contractile nanomachines that bacteria use to kill their neighbors and compete within their micro-ecosystems. PTEs help shape microbial community composition. Most PTE detection tools only detect a single PTE class. Moreover, most tailocin detection methods are largely restricted to Pseudomonas, leaving a key part of tailocin diversity uncharacterized. In this work, we present PhageTAILor (https://github.com/hjcho-bio/PhageTAILor), an integrative and fully automated pipeline that detects and classifies prophages and 3 PTE classes from bacterial genomes. PhageTAILor combines a 6-detector homology-based candidate search (geNomad, tail-gene, PHROGs-tail, SecReT6, eCIStem, and a divergence-tolerant tail-HMM detector) with a LightGBM classifier comprising 1 multiclass and 3 binary heads, trained on 6,501 bacterial genomes carrying 13,082 prophages and PTEs. A phylogeny-free feature matrix used in our model keeps predictions reproducible between model construction and user inference. PhageTAILor performs strongly at the genome level and generalizes beyond its Pseudomonas-rich training set. On a 76-strain cross-clade benchmark, PhageTAILor detected tailocins at F1 = 0.955. Furthermore, it identified 12 of 13 experimentally validated tailocins spanning five genera versus 2 of 13 for a Pseudomonas-restricted tool TattleTail. PhageTAILor also demonstrated sensitivity equivalent to viral detection tool geNomad while avoiding its higher false-positive rate. Applied to 7,925 plant- and soil-associated bacterial isolates, PhageTAILor showed that prophages in the phyllosphere and tailocins in plant-associated bacteria, whereas eCIS are enriched in soil. PhageTAILor is distributed as an open-source, modular pipeline with a command-line interface.

Protists play important roles in food chains and symbioses in soil and aquatic environments, displaying an enormous morphological and functional diversity. While most commonly found protist species are well known to science, our global-scale environmental DNA survey across soil, water, and sediments reveals dozens of novel, phylum-level phylogenetic lineages that remain to be characterized for basic morphology and function. A vast majority of these undescribed taxa occur in marine water and sediments, but some are common in soil. Most of these novel taxa have distinct substrate and habitat preferences and biogeographic patterns. To accord these lineages scientific agency and enable unambiguous scientific communication, we propose formal names for 150 species to phylum-level taxa from 25 deep lineages based on eDNA and rRNA gene long-read sequence information.

Microbial symbioses drive the evolutionary and functional diversification of eukaryotic clades, from single-celled protists to large invertebrates. However, our knowledge of host-associated assemblages (the "holobiont") is limited in microscopic animal phyla with a body size <1mm, due to practical challenges such as low biomass and difficult taxonomy of host species. Marine nematodes represent an ideal case study for rapidly advancing our knowledge of bacterial-animal symbioses, representing a globally abundant invertebrate group with strong links to terrestrial and model organism species within the same phylum. Here, we sequenced the holobionts of 220 marine nematodes and generated 815 metagenome-assembled genomes (MAGs) of host-associated bacteria/archaea. Our data indicates that 20-34% of marine nematodes harbor an obligate intracellular symbiont, often with multiple endosymbionts co-occurring within the same host. Three bacterial phyla (Pseudomonadota Bacteroidota, and Verrucomicrobiota) account for three-quarters of all nematode-associated MAGs, and the majority of these holobiont MAGs represent deeply divergent lineages in the prokaryotic tree of life. The Flavobacteriaceae (a core microbiome taxon in C. elegans and other terrestrial nematodes), were consistently recovered across phylogenetically diverse marine nematode lineages, suggesting evolutionary conservation of holobiont taxa across marine and terrestrial environments. We also report a novel chemoautotroph family (Ca. Thionematobacter) recovered from nematode hosts in both deep-sea and shallow-water habitats, and report the first confirmed instance of Cardinium endosymbionts from marine invertebrates. Finally, [~]65% of nematode-associated MAGs are able to degrade chitin, via hexosaminidase, implying that benthic invertebrate holobionts make significant contributions to global carbon cycling. These results underline the importance of evaluating symbiosis in microscopic marine invertebrates, and accelerating our understanding of animal evolution and ecosystem dynamics in vast benthic habitats.

The murine gut microbiota provides robust colonization resistance against Clostridioides difficile infection (CDI), yet murine-associated microbes remain notoriously difficult to cultivate in vitro, limiting mechanistic investigation. To identify the ecological and nutritional basis of this cultivation barrier, we leveraged CDI susceptibility as a functional readout of microbial community metabolism to infer in vivo nutrient utilization. Germ-free C57BL/6 mice colonized with varying dilutions of ethanol-treated murine microbiota were challenged with C. difficile resulting in a spectrum of CDI outcomes. Comparative metabolomics of pre-challenge fecal samples revealed a consistent carbohydrate signature: glucose accumulated in communities that resisted C. difficile challenge, whereas complex carbohydrates, including raffinose, sucrose, trehalose, lactose, sorbitol, and mannitol, were significantly depleted. The broad depletion of these complex carbohydrates supports their functional importance within the collective microbial community. Conventional glucose-based media (CMA, BHI+I, RCMT) failed to support robust growth or subculture of murine gut microbiota. Guided by the metabolomics findings, we developed Peptone Yeast Extract with Six Salts and Sugars (PYE6S), a glucose-free medium supplemented with the complex carbohydrates identified as depleted. PYE6S enabled cultivation of 22 unique Firmicutes ASVs, 82% of which lacked named cultured representatives in reference databases. These findings suggest a plausible explanation for why conventional media fail and support a metabolomics-guided framework for rational cultivation of host-associated microbiota across diverse systems. This strategy may be extended to guide media design for other host-associated microbiotas.

Determining which members of a microbial community are metabolically active remains a central challenge in microbial ecology. Although the 16S rRNA gene is the dominant marker for bacterial community profiling, it cannot reliably distinguish active cells from dormant or dead populations. As a result, complementary phylogenetic markers whose transcript abundance more closely reflects cellular activity are needed. Here, we systematically evaluated 80 Bacterial protein-coding marker genes and identified rpoB, encoding the beta subunit of bacterial RNA polymerase, as the optimal candidate. We designed a new primer pair (1528F 2041R) from a curated database of 305,274 unique rpoB sequences and validated it for quantitative PCR and amplicon sequencing of DNA and RNA templates. The rpoB qPCR assay achieved a limit of quantification two orders of magnitude lower than the benchmark 16S rRNA assay, for which a limit of detection could not be determined because of no-template-control amplification. In soil and sediment communities, rpoB recovered community composition comparable to 16S rRNA while providing a quantitative activity signal: rpoB cDNA:DNA ratios correlated significantly with taxon-level transcript abundance (R squared between 0.22 and 0.29, p < 0.001), whereas 16S rRNA ratios did not (p > 0.5). In a biological activated carbon biofilter experiment, rpoB transcript abundance tracked the decline in dissolved organic carbon removal rates across a 72 hour time series (correlation coefficients between 0.84 and 0.99), whereas 16S rRNA transcripts were uninformative (correlation coefficients between -0.4 and 0.98). These results establish rpoB as a quantitatively robust, activity-responsive complement to 16S rRNA for linking community composition to ecosystem processes.

Marine planktonic viruses play critical roles in shaping microbial communities and driving global biogeochemical cycles. However, quantitative, microbiome-wide analyses of marine prokaryoplankton-virus interactions in situ and their ecological impacts remain challenging due to the vast diversity of viral genomes and interaction modes and the limitations of existing methodologies. Here, we utilized GORG-Tropics, a global collection of 12,715 single amplified genomes (SAGs) generated from randomly sampled marine prokaryoplankton cells, to determine the frequency and modes of their interactions with viruses in the tropical surface ocean. We found 4.2% (range 1%-19% among samples) of GORG-Tropics SAGs to contain phage genomic material, with the highest frequency found in productive ocean regions. Prokaryoplankton lineages known to have high metabolic rates, including Prochlorococcus and Rhodobacterales, had a substantially larger fraction of cells associated with viruses (10-12% of SAGs) as compared to the less active but highly abundant lineages including Pelagibacterales (2.2% of SAGs). Cell-virus associations indicative of lysogeny were elevated in Alphaproteobacteria relative to other taxa. The collection of phages recovered from individual SAGs exhibited genomic diversity that bridged order-level taxonomies, indicating high diversity and genomic connectivity within wild phage populations. A substantial fraction of the observed cell-virus associations disagreed with the computationally predicted host identity of the virus, indicative of non-infective interactions. The extent of genetic exchange across tailed bacteriophages infecting different hosts, and connected to taxonomically distant phages provided further evidence for the role of non-infective phage entry in the lateral transfer between tailed phages. This study provides large-scale quantitative evidence of viral infection rates in the collective prokaryoplankton community across the global surface tropical ocean through large-scale identification and quantification of the specific phages, hosts, and modes of interaction at the resolution of individual cells. Our results confirm prior reports on the overall frequency of prokaryoplankton infections with viruses in the oligotrophic tropical surface ocean. Our findings uncover non-infective phage-cell associations that may be contributing to the lateral transfer of viral genes and the nutrition of marine prokaryoplankton.

Nitrogen fixation performed by rhizosphere bacteria has the potential to improve the sustainability of cereal crop cultivation. Deciphering the role of interspecies interactions on nitrogen fixation is crucial for devising strategies to enhance this process. To unravel the contributions of interspecies interactions, we constructed synthetic microbial communities from the bottom-up that contain diazotrophic bacteria that fix nitrogen and maize rhizosphere bacteria that do not have this capability. Interactions that impacted nitrogenase activity via growth-independent mechanisms were prevalent in the system. Nitrogenase activity increased and eventually saturated as a function of the number of inoculated diazotrophs. Using a tailored machine learning model for microbiome dynamics and explainable artificial intelligence, we deciphered species contributions on nitrogenase activity and diazotroph growth. We identified a community containing Klebsiella variicola, Herbaspirillum seropedicae, and Stutzerimonas stutzeri as a starting point for developing microbial inoculants for cereal crops. Taken together, these results provide insights into the role of interspecies interactions on nitrogenase activity.

Winogradsky columns are a classic model ecosystem for studying microbial biogeochemistry across steep gradients of oxygen and sulfide. They also remain widely used in microbiology education, introducing generations of students to microbial diversity. Yet, the genomic potential of their microbial communities remains uncharacterized. Here, we applied shotgun metagenomic sequencing to a Winogradsky column community at multiple depths, yielding 20 metagenome-assembled genomes (MAGs) representing diverse, largely uncultivated taxa. Genome-resolved analyses revealed metabolically diverse oxygenic and anoxygenic phototrophs that could potentially contribute to carbon and nitrogen fixation across all layers of the column. Most of these phototrophs also encoded one or more pathways for sulfur oxidation, which we speculated may support both energy conservation and/or sulfide detoxification by these populations. Complex carbon degradation capacity was also widespread across the MAGs, suggestive of the potential for the transformation of the column's amended organic matter (shredded coffee filters) into smaller depolymerization products and, through fermentation, organic acids. Together, these findings reveal how distinct microbial guilds might partition interconnected carbon, sulfur, and nitrogen transformations within redox-stratified systems.

Viral infection is a major yet poorly quantified driver of microbial interactions and biogeochemical fluxes in the ocean. In diatoms, which are key contributors to marine primary production, the extent to which viruses reprogram host cell metabolism and alter elemental cycling remains largely unresolved. Here, we investigated how infection by a lytic single-stranded RNA virus reshapes carbon (C) and nitrogen (N) fluxes in the ecologically relevant nanoplanktonic diatom Mediolabrus comicus. Using controlled infection experiments coupled with flow cytometry, electron microscopy, PAM fluorimetry, and stable isotope probing, we resolved infection-driven changes from the population to the cellular scale. Infection induced rapid optical shifts and cellular reorganization, including the formation of membrane-bound viral replication compartments. These changes coincided with early impairment of plastidial functions, as shown by disruption of photosystem II functionality and a concomitant decline in photosynthetic carbon fixation. In contrast, nitrogen uptake was maintained and strongly enhanced during late stages of infection, indicating sustained resource acquisition to support viral replication. This decoupling led to dynamic changes in cellular stoichiometry and, overall, to substantial reductions in population-level carbon and nitrogen assimilation due to growth inhibition. Together, these findings demonstrate that diatom RNA virus infection reshapes host carbon and nitrogen metabolism, with cascading effects on elemental cycling. Our results identify diatom RNA viruses as important drivers of marine biogeochemical processes, with implications for primary production and the fate of organic matter in the ocean.

The SAR11 clade, also known as the order Candidatus Pelagibacterales, is among the most abundant bacterial lineages in the ocean and plays central roles in marine biogeochemical cycles. However, many SAR11 genes remain functionally uncharacterized, highlighting the need for a comprehensive, integrated catalog that supports genomic, functional, and ecological analyses across the clade. Here, we present the SAR11 Genome Atlas, an interactive ortholog group (OG)-centered web resource that integrates 542 SAR11 genomes, including all 132 cultured strain genomes, with functional annotations, synteny, phylogenetic distribution, metatranscriptomic expression, and predicted protein structure information. To demonstrate its utility, we used environmental expression profiles to identify OGs associated with high-latitude environments, recovering OGs known to be involved in cold adaptation and proposing a hypothesis for the function of uncharacterized protein. We further analyzed phylogenetic distribution patterns to identify mutually exclusive functional modules, including candidate alternative systems for Mn/Zn homeostasis and phosphate acquisition, and to associate these modules with distinct oceanographic environments. Together, these case studies demonstrate that the SAR11 Genome Atlas supports complementary analyses that connect environmental signals to genes of interest and use phylogenetic or functional distributions to generate hypotheses about ecological specialization. Through a user-friendly web interface, the SAR11 Genome Atlas enables researchers to explore genomic, environmental, and structural information without specialized computational expertise. All data and analysis outputs are freely accessible online at [https://stsnsn.github.io/SAR11_Atlas/]. The SAR11 Genome Atlas thus provides a scalable framework for generating and testing hypotheses that connect SAR11 genomic variation to protein function and oceanographic context, supporting advances in marine microbial ecology and biogeochemistry.

Overproduction of indole-3-acetic acid (IAA) by rhizosphere bacteria disrupts plant auxin homeostasis and induces root growth inhibition (RGI). Variovorax reverses this effect by degrading IAA, but the underlying pathway remains incompletely resolved. Here, using genetics, metabolomics, and isotope tracing, we define a nine-gene region (iadCDEFGHIJK2) required for IAA catabolism in Variovorax paradoxus CL014, identify a previously uncharacterized intermediate (CHNO), and revise the early oxidative steps of the pathway. Although IadDE adopts a Rieske-type oxygenase architecture, our data are consistent with the IadCDE complex functioning as a monooxygenase during IAA degradation. Using these insights, we introduce iad genes into two rhizobacterial chassis, Polaromonas MF047 and Paraburkholderia MF376. These engineered strains degrade IAA and alleviate RGI induced by exogenous IAA, an auxin-producing strain, and a synthetic bacterial community. Engineered Paraburkholderia MF376 delivers the strongest performance, improving plant growth in natural soil. Together, these results establish a framework for engineering auxin-balancing root commensals.

The Southern Ocean harbors diverse marine microbial communities shaped by both local oceanographic conditions and dispersal limitations. However, this knowledge is mainly based on coastal Antarctic sites, whereas circumpolar Antarctic open sea and subantarctic ecosystems remain poorly explored. Here, we characterize marine microbial communities (using 16S rDNA high-throughput sequencing) and marine oceanographic data across two regions: the Subantarctic, involving two localities (the Magellan Strait and the Beagle Channel), and Antarctic open sea, involving two localities (Eastern Indian and Central South Pacific). We found extensive differences across regions and localities, characterized by distinct taxonomic patterns, alpha diversity, microbial composition, and enriched taxa profiles. Despite these differences, Clade Ia, Amylibacter, NS5 marine group, and NS2b marine group exhibited high prevalence across regions. Oceanographic parameters had variable relationships with microbial alpha diversity across regions: Sea surface temperature and salinity had a negative and positive correlation, respectively, in the Magellan Strait during 2024. In the Antarctic region, dissolved oxygen displayed a negative correlation in the Indian Ocean during 2024, whereas salinity displayed a more variable relationship in the Indian Ocean: positively correlated during 2024, while negatively correlated during 2025. Collectively, our results highlight a strong microbiological biogeographic structure in the Southern Ocean, both across broad scales (between Subantarctic and Antarctic regions) and within regions. Furthermore, our results show dynamic relationships between oceanographic variables and marine microbial diversity across Antarctic and Subantarctic regions.

Carboxylases are among the most important enzymes in nature as they catalyze the fixation of inorganic carbon (CO2), a central step in the global carbon cycle. In addition to their well-known function in autotrophic CO2 fixation, many carboxylases play a role in the heterotrophic assimilation of organic compounds. In this study, we provide genomic evidence for an assimilatory carboxylation pathway involved in acetone degradation in photoheterotrophic bacteria from metagenomes collected along a latitudinal transect of the Arctic Ocean. This curious metabolism was linked to a single population of Gammaproteobacteria (Porticoccus arcticus). P. arcticus has a streamlined genome compared to Porticoccus relatives but has maintained a complete acetone carboxylation pathway while acquiring multiple proteorhodopsin genes by lateral gene transfer. Arctic Ocean metatranscriptomes revealed the acetone carboxylase and rhodopsins genes were among the most highly expressed P. arcticus genes in oligotrophic Arctic surface waters. P. arcticus sequences were consistently detected, and often abundant (up to 9%), in a multiyear Arctic Ocean 16S rRNA time-series, supporting its ecological significance in Arctic marine systems. Overall, this work reports a metabolic module (acetone carboxylation) in the ocean that may allow photoheterotrophic bacteria to enhance their biosynthetic capacity via CO2 assimilation.

While the suppression of toxic cyanobacteria by aquatic plants has long been recognized, few studies have clearly differentiated between the allelopathic effects of the plant itself and the inhibitory influence of its associated microbiome. This study aimed to clarify the primary inhibitory agent by pre-culturing Myriophyllum spicatum (Eurasian watermilfoil) under three conditions: (1) BG11 medium, (2) live Microcystis aeruginosa KW culture, and (3) a Microcystis-symbiotic microbiome (excluding Microcystis cells). After a 7-day pre-culture, Myriophyllum shoots were transferred to fresh Microcystis culture. The Myriophyllum pre-cultured in Microcystis culture exhibited rapid inhibition against Microcystis (84% within day 1), whereas the Myriophyllum pre-cultured in BG11 medium showed delayed responses (89% by day 7). In contrast, inhibition remained below 50% in the Myriophyllum pre-cultured with the Microcystis-symbiotic microbiome. Notably, plant-derived soluble compounds exhibited weak inhibitory effects, whereas the microbiome showed stronger inhibitory activity, indicating that the plant-associated microbiome plays a more dominant role than the plant itself. Exposure to Microcystis triggered significant shifts in plant-symbiotic microbial community composition, leading to rapid enhancement of inhibitory activity in the Myriophyllum microbiome. Microbial community analysis identified 28 bacterial taxa closely associated with the inhibitory response, including strains involved in organic matter degradation, adhesion, biofilm formation, and predatory behavior. Meta-transcriptomic analysis further confirmed increased expression of genes related to bacterial adhesion, biofilm formation, and carbohydrate metabolism following Microcystis exposure, highlighting functional adaptations linked to cyanobacterial suppression. These findings underline the role of microbiome-mediated cyanobactericidal mechanisms, providing new insights into a nature-based solution for mitigating Microcystis-dominated harmful algal blooms.

Anaerobic gut fungi (AGF) are central to the degradation of plant material in the digestive systems of herbivores. However, how their environment influences their colonisation and degradation of complex biomass is unclear. Here, cellulose filter paper was used as a simplified model of the plant cell wall to investigate how the presence of free sugars in the rumen can affect AGF growth and degradative responses of phylogenetically distinct AGF isolates. From this, galactose was revealed to be inhibitory to both Neocallimastix frontalis and Caecomyces communis, and mannose inhibitory to C. communis. Complete inhibition of C. communis growth was conserved when galactose and mannose were added in their polymeric forms, whereas in contrast, N. frontalis growth was unaffected. This indicates, depending on the AGF isolate, the presence of free sugars and their polymeric form may influence AGF growth through regulatory and metabolic interactions - even if the sugar cannot be utilised for growth as the sole substrate. Collectively, this work highlights the functional diversity in AGF carbohydrate responses and the need for greater understanding of their metabolic regulation for applications in lignocellulosic bioconversion and ruminant nutrition.

Giant viruses (GV) are increasingly recognized as important ecosystem regulators. While metagenomics has uncovered extensive GV diversity, the global distributions of individual species and the biogeographic processes driving the pattern remain poorly understood. Here, we reconstructed GV metagenome-assembled genomes (MAGs) from 35 globally distributed deep freshwater lakes spanning five continents, aiming to identify their biogeographic patterns. The resulting 1663 non-redundant MAGs significantly expanded the known freshwater GV diversity, with [~]84% lacking a previously reported species representative. These MAGs were grouped into cosmopolitan and geographically restricted lineages. We identified 27 cosmopolitan GV species spanning multiple viral lineages, including families of Imitervirales, Pimascovirales, and mirusviruses order Styxvirales. The cosmopolitan species were characterized by their larger genomes and expanded gene repertoires of host-interaction functions, which may facilitate interactions with diverse hosts and contribute to their global distributions. The presence of geographically restricted species and the stronger distance-decay in community similarity observed in freshwater than marine ecosystems suggest that physical connectivity between ecosystems is an important factor influencing GV dispersal. We identified 312 and 177 GV MAGs almost exclusively associated with the epilimnion and hypolimnion, respectively. This water-layer preference of individual MAGs was highly consistent across lakes, suggesting conserved vertical partitioning in association with the thermal stratification of the water column. Overall, our findings reveal that GV biogeography in deep freshwater lakes is structured by the combined influence of horizontal dispersal limitation, vertical partitioning, and lineage-specific evolutionary histories.

Predatory protists are important in shaping terrestrial microbial ecosystems, but their roles in the phyllosphere, or the communities on aerial plant surfaces, are poorly understood. Previous work found that the order Colpodida dominated heterotrophic protist communities in the phyllosphere. While most protists were sporadically present, a few Colpodida variants were prevalent and abundant, indicating that these variants may represent species adapted to the phyllosphere. To identify these organisms, we cultured colpodids from field-collected tomato leaves and performed phylogenetic analysis of the 18S rRNA gene. Five of nine independent isolates matched the most prevalent Colpodida variant previously identified as leaf-enriched through amplicon sequencing, and these isolates comprised a novel clade of Paracolpoda steinii. When compared to a maize root isolate of Colpoda inflata, an abundant rhizosphere ciliate, a P. steinii isolate was similar in size and growth yield on E. coli, but grew to higher yields and formed large cyst clusters when incubated with model phyllosphere bacteria prey Erwinia and Pseudomonas. We developed and validated quantitative PCR (qPCR) methods for detection and cell abundance estimation of the P. steinii phyllosphere clade, C. inflata, and the order Colpodida in environmental samples. In inoculated greenhouse plants, qPCR-estimated protist populations matched measured inoculum levels, and protist inoculum was still detectable after five days. In an uninoculated tomato field, P. steinii was detected on all plants, with greatest abundances observed in lower leaves and after a rain event. P. steinii comprised up to 18.7% of total leaf Colpodida populations, which were estimated at up to [~]1400 organisms per gram of fresh weight. The findings demonstrate that Colpodida communities are consistently present on tomato leaves, dynamically affected by the abiotic environment, and include significant populations of P. steinii. We propose that the P. steinii isolates and qPCR tools presented can be used as a model system to investigate colonization and distribution patterns, biotic interactions, genetic adaptations, and agricultural applications of leaf predation.